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We show how to use muhiple tracers of large-scale density with different biases to measure the 
redshift-space distortion parameter {3 = f = b~^dlnD/d\na (where D is the growth rate and a 
. . . the expansion factor), to a much better precision than one could achieve with a single tracer, to an 

arbitrary precision in the low noise limit. In combination with the power spectrum of the tracers 
t*—^ ' this allows a much more precise measurement of the bias-free velocity divergence power spectrum, 

' f^Pm - in fact, in the low noise limit f^Pm can be measured as well as would be possible if velocity 

divergence was observed directly, with rms improvement factor ~ [5.2 (/^^ +2/3 + 2) /P^] ^''^ (e.g. 
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10 times better than a single tracer for fi = 0.4). This would allow a high precision determination 
of fD as a function of redshift with an error as low as 0.1%. We find up to two orders of magnitude 
improvement in Figure of Merit for the Dark Energy equation of state relative to Stage II, a factor 
' of several better than other proposed Stage IV Dark Energy surveys. The ratio 62/&1 will be 

determined with an even greater precision than (3, producing, when measured as a function of 
scale, an exquisitely sensitive probe of the onset of non-linear bias. We also extend in more detail 
1-^ ' previous work on the use of the same technique to measure non-Gaussianity. Currently planned 

redshift surveys are typically designed with signal to noise of unity on scales of interest, and are 
Q , not optimized for this technique. Our results suggest that this strategy may need to be revisited as 

■ there are large gains to be achieved from surveys with higher number densities of galaxies. 

^ _ PACS numbers: 98.65.Dx, 95.35.-fd, 98.80.Es, 98.80.-k 

> . 

, I. INTRODUCTION 

(N 

' Growth of structure in the universe has long been recognized as one of the most powerful ways to learn about the 
nature of dark energy and other properties of our universe. Currently the most promising method is weak lensing 
tomography, which traces the dark matter directly and can measure the growth of structure by sphtting the source 
galaxies by their (photometric) redshift. For example, the Dark Energy Task Force concludes that of all proposed 
, next generation experiments weak lensing holds the best promise to succeed in improving our knowledge of dark 
energy However, weak lensing is not necessarily the ideal method: it measures the structure projected over a 
broad window in redshift and its abihty to probe rapid changes in growth rate is limited. Moreover, one is measuring 
• a 2-dimensional projection only as opposed to the full 3-dimensionaI information, drastically reducing the amount of 
J-j ' available information. In addition, for all of the methods proposed so far the requirements on systematic error exceed 
. . 1 what is achievable today, so it is worth pursuing multiple methods until we have a better control of systematics. 

Galaxy clustering has been the favorite method of measuring large scale structure in the universe and it is likely 
this will continue also in the future. The main reason is that galaxies are easily observed and that by measuring 
their redshift one can reconstruct the 3-dimensional clustering information, in contrast to weak lensing or cosmic 
microwave background anisotropics which only measure a 2-dimensional projection. The relation between galaxy and 
dark matter clustering is however not straight-forward. In the simplest model of linear bias galaxies trace dark matter 
up to an overall constant called linear bias. The bias cannot be predicted from the theory and as a result galaxy 
clustering alone cannot measure the growth of structure with redshift. If one could determine bias with sufficiently 
small error then galaxy clustering would become the leading method due to its higher information content. 

There are several methods proposed to determine the bias. One is to use weak lensing, specifically galaxy-galaxy 
lensing which measures the cross-correlation between galaxies and dark matter. This is proportional to bias b and in 
combination with the galaxy auto-correlation function which scales as one can eliminate the dependence on bias 
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to measure fi^^pPm (fc) , where ^Imo is the matter density parameter today and Pm{k) is the matter power spectrum 
at a given redshift [2, [Sj . Alternatively, one can also measure the halo mass distribution with weak lensing which, 
in connection with the theoretical bias predictions, can determine the bias and thus Pm{k) A second method to 
determine the bias and thus Pm{k) is to measure the three-point function A third method, and the one we focus 
here, is to measure the redshift space distortion parameter f3 = b~^f = b^^d\nD/d\na (where b is the bias of the 
galaxies, D is the linear theory growth factor, and a is the expansion factor) 0,0, BUI. In combination with the galaxy 
power spectrum this gives f^Pm- Note that these methods give somewhat different dynamical measurements once 
the bias is eliminated, so to some extent they are complementary to each other. However, none of these methods is 
presently competitive in terms of derived cosmological constraints, as they all have rather large statistical errors from 
the current data, although this may change in the future as data improve and new analysis methods are developed 




In this paper we focus on redshift space distortion parameter /3 as a way to determine the bias. Clustering of 
galaxies along the line of sight is enhanced relative to the transverse direction due to peculiar motions and this allows 
one to determine (3. Current methods require one to compare the clustering strength as a function of angle relative 
to the line of sight, but this method is only applicable on large scales where linear theory holds. As a result, the 
sampling variance limits its statistical precision. Recently, 14] proposed a new method for measuring the primordial 
non-Gaussianity parameter /nl, by comparing two sets of galaxies with a different bias, and a different sensitivity to 
/nl [IBIj which allows one to eliminate the sample variance. Optimization of this technique was investigated by p^ . 
Here we apply this technique to the measurement of the redshift-space distortion parameter, which in turn improves 
the measurement of the velocity divergence power spectrum Pgg = pPm{k). As we show here, this approach can 
in principle measure /3 perfectly and Pgg as well as if we observed velocity divergence directly. If this promise were 
realized from the data it would allow for a much higher statistical power than weak lensing or other 2-dimensional 
projections. We also show that our method, and large-scale structure surveys in general, become even more powerful 
when additional cosmology dependence, such as the Alcock-Paczynski effect is included. 

We begin by presenting the basic method, followed by an analysis of the expected improvement as a function of 
survey parameters. This is followed by predictions for some of the existing and future surveys in terms of expected 
improvement of dark energy parameters. On small scales the deterministic linear bias model eventually becomes 
inaccurate and for this reason we often quote results as a function of the maximum usable wavenumber, fcmax. While 
we have some idea what its value should be, we leave a more detailed analysis with numerical simulations for the 
future. One should note, however, that the multiple-tracer approach is likely to actually help in disentangling non- 
linear bias effects on quasi-linear scales, making those scales potentially much more useful than they would otherwise 
be. We conclude with a summary and a discussion of future directions. 



II. METHOD AND RESULTS 

The density perturbation, 5gi, for a type of galaxy i, in the linear regime, in redshift space, is [Tsj 

Sg^ - + /m') 5 + (1) 

where bi is the galaxy bias, /i — 5 is the mass density perturbation, and is a white noise variable which 

can represent either the standard shot-noise or other stochasticity. All equations are understood to apply to the real 
or imaginary part of a single Fourier mode unless otherwise indicated. In this paper we will consider two types of 
galaxies, type 1 with bias 6, and type 2 with bias ab. The perturbation equations can then be written 

5gl = f{P~^+^l^)6 + e^, (2) 

and 

5g2^f{a^^+^l^)S + e2. (3) 

We are denoting (} = f /b (the equivalent distortion parameter for the 2nd type of galaxy is (5/ a). 

The traditional method to determine (3 is to look at the angular dependence of the two-point correlation function 
or its Fourier transform, the power spectrum. The correlations will be enhanced along the line of sight (where /i = 1) 
relative to the direction perpendicular to it (where ^ — Q), but to observe this enhancement one must average over 
many independent modes to beat down the sampling (or cosmic) variance. This is because each mode is a random 
realization of a Gaussian field and there will be fluctuations in the measured power even in the absence of noise. By 
combining a measurement of (3 with that of the galaxy power spectrum one can determine f^Pm, which no longer 
depends on the unknown bias of the galaxies. This method has been applied to the data, most recently in ^, Ji, i8|, i9| , 
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and is limited by the accuracy with which we can determine (3. For example, for the analysis in which currently 
has the highest signal to noise measurement of (3, the error on the overall amplitude of galaxy power spectrum Pgg 
is about 1% adding up all the modes up to fc = 0.1 /iMpc^^, while the error on (3 is about 12%, so the error on 
reconstructed pPm = P^Pgg is entirely dominated by the error on f3. Predictions for the future surveys with this and 
related methods can be found in [H, [l3| ■ 

To understand the main point of the new method we are proposing here let us consider the situation without the 
noise, which would apply if, for example, we have a very high density of the two tracers sampling the field, and no 
stochasticity. In that case we can divide equation [3] by equation [3 above to obtain 



Sg2 
Sgl 



a/3- 



(4) 



This expression has a specific angular dependence, allowing one to extract a and (3 separately. Note that there is no 
dependence on the density field 5. The random nature of the density field is therefore not affecting this method and 
we can determine (3 exactly in the absence of noise. More generally, the precision with which we can determine (3 is 
controlled by the (shot) noise, i.e. density of tracers, rather than the sampling variance. In order to address the gains 
in a realistic case we must perform the full analysis, which we turn to next. 

Generally, the noise variables can be correlated with each other, although they would not be for standard shot-noise. 
For example, [l9l | showed that one generically expects non-linear structure formation to generate a full covariance 
matrix for the noise at some level, so this will probably be the ultimate limit for this kind of measurement. The 
covariance matrix of the perturbations is 



C . 



{5gl&g2) 



Pee 
2 



N 



(5) 



where Peg = 2p 1^5'^^ and iVy = 2 {ei^j) (note that Pee and N are the usual power spectrum and noise - the factor 
of 2 comes from the fact that 5 and e are only the real or imaginary part of a Fourier mode). If we assume that 
the noise matrix is known then the covariance matrix for two types of galaxies is a function of three parameters: 
the velocity divergence power spectrum amplitude. Pee, the redshift-space distortion parameter, (3, and the ratio of 
biases a. We work with a instead of a second (3 parameter because the ratio of biases will generally be determined 
substantially more precisely than /3, which means that the measurement of the second 13 parameter would be almost 
perfectly correlated with the first. The Fisher matrix for the measurement of these parameters is 



^AA' = ^Tr [CaC-^CvC-i] 



(6) 



where Ca = dC/dX and A are the parameters. 

For any single mode, the inverse of the Fisher matrix is singular, i.e., we can only constrain two parameters, not 
three. However, adding Fisher matrices for modes with different /i breaks this degeneracy. Generally, the total Fisher 
matrix will be an integral over modes with all angles. As a simple example, we assume first a pair of modes with 
/i = and /i = 1 and compute the error on the parameters in the small noise limit: 



— ^ — Xii — 2Xi2 



22 



(7) 



where — NijlbibjP„ 



1£ 



a^l + f]f + {a + (3f 



Xn-2 



(1 +py+a{l + 13) (a + [3) + 2c? (1 + i3y X; 



f3^ {a-iy 



(8) 



and 



(9) 



The key points are that the only lower limit on the errors on a and /? from this single pair of modes is set by the 
achievable noise-to-signal ratios on the tracers, Xij, and the error on Pge is the error one would obtain for a simple 
power spectrum measurement from two modes, with no degradation due to degeneracy with the bias parameters (this 
is of course only true to leading order in the small noise limit). 
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In contrast, for a single type of galaxy, for the same pair of radial and transverse modes, in the low noise limit we 
, J. 

(10) 



find 

^1 _ (1 + /3) 
and 



< _ 2(/32 + 2/3 + 2) 

i.e., here the signal-to- noise per mode on /3 saturates at greater than unity for low noise, and the overall power 
measurement is substantially degraded because we cannot determine (3 very accurately and this then limits the 
precision with which the velocity divergence power spectrum can be determined as well. If we assume the first tracer 
is much better sampled than the second, we get the ratio of errors 



cr| (2 tracers) ^ 2a^X22 _ 2N22 
aj (1 tracer) ~ {a - if " {b2 ~ b^ f P„ 

and 

cr|,^^ (2 tracers) f}"^ 



al (1 tracer) 2 (/32 + 2^ 2) 



(12) 



(13) 



We see that the improvement in the measurement of (3 can be arbitrarily large in the limit N22 — > 0, while the 
improvement in Pgg is limited by sample variance but is generally quite large, e.g., a factor 10^/^ (rms) for /3 = 1, and 
50^/2 for f3 = 1/3. All of these formulae are only good as long as the noise is small relative to the power spectrum of 
the difference between the two tracer fields (obviously the two tracer case should never be worse than the one tracer 
case). 

A yet more realistic analysis must integrate over all modes rather than look just at the two specific modes as done 
above. In doing so we find that the two-mode approximation above substantially underestimates the improvement 
in the Pee measurement. If we integrate the Fisher matrix over the full range of /j,, the improvement is a factor of 
approximately 2.6^/^ larger than above, e.g., a factor of 25^^^ for (3 — 1, or 138^/^ for /3 = 1/3. The result of the ji 
integration in the idealized case of two perfect tracers remains simply: 

^ = - (14) 
Pee' ^'^> 

where N is the effective number of discrete modes sampled (counting real and imaginary parts separately). Thus, the 
combination of two perfect tracers produces the measurement of /^P^ that one would obtain if one could observe the 
velocity divergence field directly and perfectly. 



A. General usefulness of redshift-space distortion measurements 



Of course, we will never have perfect sampling, so it is necessary to numerically evaluate the exact Fisher matrix 
with finite noise to determine the improvement for real surveys. We assume the survey volume, V, is roughly spherical, 
so the Fisher matrix for a survey will be 

Fv-^J^ "^^'^ k^dk J d^F {k, fi) (15) 

where F (fc, fi) is computed for a single mode using Eqs. (O and ([6]) (remembering that both real and imaginary parts 
must be added). We assume the minimum usable k is fcmin — 2n/V^^^. 

In Figure [1] we show the errors on the amplitude of the velocity divergence power spectrum Pgg that one can 
obtain from a generic 3/4-sky (30000 square degree) survey, as a function of redshift. This may seem unreasonably 
optimistic/futuristic, but it may be possible surprisingly quickly if the 21 cm emission intensity mapping method is 
successful [2^, as discussed further below. We assume either k^a^{z) = 0.1 [D (z) /D{0)]~^ hMpc^^, or kjnax{z) = 
0.05 [D (z) /D (0)]"^ /iMpc"\ with the redshift dependence motivated by [Uill. We give results for different values 
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FIG. 1: Projected fractional error on the normalization of Peg = f^Pm, for 30000 square degrees, in redshift bins with width 
dz = 0.2. The upper and lower green (short-dashed) lines show the constraints from single tracers with fo = 2 and h = 1, 
respectively. Black (solid) lines show the two tracers together, each with, from top to bottom, S/N=l, 3, 10, 30, 100 at 
k = 0.4/iMpc~^. The red (long-dashed) line shows the case where both tracers are perfectly sampled. For the left panel we 
assume fcmax(2) = 0.1 [D {z) /D (0)]"^ /iMpc"\ while the right assumes /cmax(2) = 0.05 [D {z) /D (0)]"^ /iMpc"\ 



of Xij ai k = (for i ~ j, i.e., the off-diagonal element is zero). Surveys with S/N~ 1 at fc = 0.4ft,Mpc~^ 

recover ~ 90% of the total possible information from BAO measurements, i.e., the S/N ratios we give should be 
interpreted as being relative to a nearly ideal BAO survey. In these calculations we are assuming the shape of the 
power spectrum is known so that we can compress the information from all the modes into a single number, that of 
the power spectrum amplitude. This is discussed further below. 

We see from Fig. [T] that one obtains large improvements using the multi-tracer method, even for S/N as small 
as 3. For example, at z ^ 1, with fcinax(2) = 0.1 [D (z) /D{0)]~ hMpc~^, the gain is equivalent to increasing the 
survey volume by a factor of 3 over a single, perfectly sampled, unbiased tracer. For the more conservative kj^axiz), 
the results are substantially degraded, but the multi-tracer method actually becomes more valuable, at fixed absolute 
S/N level, because the S/N level in the usable range of k is higher. 

It is interesting to note that, in the single-tracer case, low bias is substantially better than high bias. This is easy 
to understand from Eq. [T] - lower bias simply means that power due to redshift-space distortions makes a larger 
fractional contribution to the total. The practical consequence of this observation is that, once BAO surveys using 
high bias objects have been completed, an opportunity exists for improvement by observing lower bias objects at 
similarly low S/N, i.e., even if the S/N is not large enough to decisively exploit the multi-tracer method in this paper. 

To put these results into a broader context. Fig. [5] shows how velocity divergence power spectrum Pgg — pPm 
depends on the underlying cosmological parameters of interest, as a function of redshift. The derivatives with respect 
to parameter p^, dlnPge/dpi, are taken at fixed values of the other parameters in the following set: cOm = ^mh^, the 
matter density; ujb = ^bh^ , the baryon density; 0^, the angle subtended by the sound horizon at the CMB decoupling 
redshift; and w', defined by w{z) = + (1 — a) = p{z)/ p{z) where here p{z) and p{z) are the dark energy 
pressure and density; fifc, the effective fraction of the critical density in curvature; and A and rig, the amplitude and 
power law slope of the primordial perturbation power spectrum. Fixing a;,„, w?,, ^s, A, and guarantees that changes 
in wq, w' , and Qk have minimal effect on the CMB. The scale for the parameter dependences is chosen using the error 
bars on each parameter that we project for the scenario described for Fig. [5l including the Planck CMB experiment, 
and Pgg and BAO constraints up to Zmax = 1- While it is easy to see that Pgg measurements have a lot of power to 
constrain parameters, especially the dark energy equation of state, exactly what constraints can be obtained is not 
obvious to the eye, because one needs to account for degeneracy between parameters. 

In the derivatives in Fig. [21 and the following calculations, we assume that the constraints on Pgg can be taken as a 
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FIG. 2: Dependence of Pee {k — 0.1 /i Mpc"'^, on cosmological parameters. The derivatives are at fixed values of the other 
parameters, including not-shown parameters ujb, ds, and Us- The derivatives are normalized by typical projected errors on the 
parameters, for the scenario described in Fig. (5] with Zmax ~ 1. Lines correspond to: black (solid): wq, green (long-dashed): 
w' , red (dot-short-dashed): fifc, magenta (dot-long-dashed): uim, cyan (dotted): log A. Error bars show the projected errors on 
Pee for k^^^z) = 0.1 [D (z) /£'(0)]"' hMpc-\ and S/N = 10. 



measurement of Pgg at fc = 0.1 hMpc~ , for all z and all cosmological models. This allows us to apply the constraints 
directly to the product f{z)D{z), with the power spectrum entering only through a single overall normalization. We 
do this for simplicity and transparency (i.e., so one can see clearly where the constraints are coming from), but it 
is imperfect in two ways: first, the measurement at different redshifts is really weighted toward different fc, so the 
shape of the power spectrum enters in converting between them; and, second, we can not measure the power spectrum 
in h^^ Mpc units, only radial velocity separations and angular separations, which leads to sensitivity to the Hubble 
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FIG. 3: Improvement in Dark Energy constraints (quantified by the DETF FoM), relative to Planck + DETF Stage II, 
when the Pee constraints from Fig. [1] are added. The lines refer to the same cases as in Fig. [TJ and again the left panel is 
fcmax(2)=0.1 [D(z)/D(0)]"^ /i Mpc" ^ while the right is fc^ax («)= 0.05 [D (z) / D (0)]''^ hMpc~\ 



parameter, H{z), and angular diameter distance, Da{z) (this includes, but is not entirely, the Alcock-Paczynski effect 
|17|). We will revisit these issues later. 

In Fig. [3] we show the contribution these measurements of Pgg can make to the study of dark energy, quantified by 
the Dark Energy Task Force (DETF) Figure of Merit (FoM) , proportional to the inverse of the area within the 95% 
confidence contours describing constraints on wq and w' [1]. We see that there is the potential for a huge improvement 
in dark energy constraints, by two orders of mag nitude in the FoM, when fcmaxC^) = 0.1 [D [z) /D (0)] /iMpc"\ 
well beyond even the optimistic Stage IV constraints envisioned by the DETF (which gave an improvement of a factor 
^ 20 in the units of this figure). The FoM improvements for the fcmax('Z) = 0.05 [D (z) / D (0)]^ /iMpc^^ constraints 
are substantially weaker of course, but still potentially spectacular. 

Finally, we add the BAO distance measurements that one would obtain from the same surveys, with the resulting 
FoM shown in Fig. [4l We see that even a minimal Peg measurement adds substantially to the BAO-only measurement. 
To give more meaning to these FoM improvement numbers: in the best case shown here, the constraint on w{zp) 
at the pivot point Zp ~ 0.8 (where the errors on w{zp) and w' are uncorrelated) is ±0.0023, with w' is measured to 
±0.031. In the less optimistic case of S/N — 10 out to Zmax = 2, with fcniax(^) = 0.1 [D (z) / D {0}] ^ ZiMpc"^, we 
find w{zp = 0.45) is constrained to ±0.0074, and w' to ±0.091. 

Finally, in Fig. [5] we show the constraints on all of the parameters where Pee helps significantly, not just dark 
energy. This figure is for the specific case with Planck, Stage II, BAO from the same galaxy survey included, 
fcmax(^) = 0.1 [D (z) /D{0)]~^ /iMpc~^, and S/N=10 (note that, with one exception that we will indicate, fcmax 
generally does not apply to the BAO calculation, which essentially follows the procedure in [13] )■ We see that 
dramatic improvements will be made in the constraints on Qk, (the dark energy density), and even and A. 



B. Specific example surveys 



We now explore a few planned or existing surveys. 
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FIG. 4: Improvement in Dark Energy constraints (quantified by the DETF FoM), relative to Planck + DETF Stage II, when 
the Pee constraints from Fig. [T]are added, and the BAO distance measurements from the same survey is added. The lines refer 
to the same cases as in Fig. [TJ except the blue (dotted) line is the case with BAO alone. The left and right panels again show 
the stronger and weaker values of femax, respectively. 



1. High redshift, high volume space mission 

As a first example, we consider scenarios motivated by the EUCLID mission [2^ (EUCLID is the combination of the 
missions formerly known as SPACE and DUNE - there is a proposal similar to SPACE in the USA named ADEPT). 

The proposal is to cover 1/2 of the sky over a wide range of redshifts. Our calculation considers 100 (/i~^Gpc) 

centered at z = 1.4. We use the planned galaxy density 0.0016 {h~^ Mpc) ^ (L. Guzzo, private communication). We 
assume a mean bias of 6 = 2 (somewhat arbitrarily ~ the bias could be significantly higher if one is able to select 
galaxies in the most massive halos at that time). Fig. [6] shows the f^Pm measurement for various modifications or 
additions to the EUCLID mission. We see (black, solid line, in the upper left) that this survey is a long way from 
achieving the best possible measurement of Pge, by almost a factor of 10 in the rms error. Viewed as a single tracer, the 
problem is not noise (compare to the magenta, dotted line), it is simply that, for 5 = 2, we lose this much constraining 
power to degeneracy with the bias. A single perfectly sampled tracer with b — 1 would do better, but still be far from 
as good as possible (blue, dot-dashed). If we hypothetically split the sample into b = 1.5, fi = 0.0008 { h^^ Mpc) 

and 6 = 2.5, n = 0.0008 {h^^ Mpc) subsamples (chosen for illustration, not because we know this is possible), in 
order to exploit the method of jtU], the gains are non- negligible, especially on large scales where they can more than 
double the effective volume of the survey, however, the noise is too large to achieve the much larger gains that are in 
principle possible. If we could leave the EUCLID-like sample intact and add a perfect unbiased tracer as the second 
field, we would get a substantial gain from the combination, equivalent to multiplying the volume of the survey by a 
factor of ~ 10 — 17 (the exact improvement factor depends on fcmax), although a factor '--^ 3 of this would be obtained 
from the perfect unbiased tracer alone. As a practical matter, a factor of ten larger number density for the unbiased 
tracer over the biased tracer is nearly equivalent to perfect, i.e., the rms error is within ^ 10% of the infinite density 
error. 

The ratios shown in the upper left of Fig. [6] are largely independent of the volume of the survey. The upper right of 
Fig. [S] shows how well one can do in an absolute sense with a 100 cubic Gpc/h survey. One critical point to take away 
from this figure is how vital it is to push fcmax to be as large as possible. The overall error on a power spectrum scales 
like fcj^ax iong as one has good S/N, so a mere 30% slippage in fcmax (e.g., O.l/iMpc^^ instead of 0.13/iMpc~^) 
is equivalent to throwing away more than half of the survey volume! Similarly, when one is using the multi-tracer 
method in this paper, increasing fcmax by a factor of two is roughly equivalent to an order of magnitude increase in 



9 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




12 3 4 

z 



FIG. 5: Improvements in constraints on various parameters, vs. «max, for the case including BAO, fcmax(z) = 
0.1 [D (z) / D {0)]~^ ^Mpc~^, and S/N = 10 for the Pee measurement. Lines are: black (solid): wo, green (long-dashed): 
w' , red (dot-short-dashed): Qk, blue (dashed): Qa, magenta (dot-long-dashed): corn, cyan (dotted): logyl. 

the number of galaxies. Future analyses will need to explore fcmax of the single versus double tracer method, but for 
this numerical simulations are needed, which are beyond the scope of this paper. 

Figure [6] also shows the projected error on f3, which shows very similar dependences to pPm- In fact, the (3 and 
f^Pm errors are almost perfectly correlated, to the point where it is probably most accurate to think of these surveys 
as making a very precise measurement of b^Pm which is converted to f^Pm by multiplying by /3^, with the resulting 
error dominated by the error on (3. Note, however, that with very low noise, the /3 error will continue to decrease 
toward zero, while the b^Pm error will not and therefore will eventually dominate the f^Pm error (of course, this is 
just the cosmic variance limit for a direct measurement of velocity divergence). Fig. [S]also shows the projected error 
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FIG. 6: Upper left panel: Projected error (inverse variance) on the normalization of Pee = f^Pm, relative to the error one could 
achieve by observing velocity divergence directly (which is also obtainable by observing two noise-free tracers with different 
bias), as a function of fcmax, for the EUCLID mission-motivated example described in the text. The black (solid) line shows 
EUCLID-like galaxies treated as a single tracer. The magenta (dotted) and blue (dot-short-dashed) lines show single perfectly 
sampled tracers with b — 2 and 6=1, respectively. Green (short-dashed) shows the EUCLID-like galaxies split into two groups 
with different bias, as described in the text. Red (long-dashed) shows the EUCLID-like galaxies plus a perfect unbiased tracer, 
while cyan (dot-long-dashed) shows a density n = 0.03 ( Mpc) plus an unbiased tracer. Where there are two lines of the 
same type, the upper one uses a genuinely perfectly sampled tracer, while the lower one shows density n = 0.03 ( Mpc) ^ . 
Because the error is always dominated by k's near fcmax, this plot (and subsequent plots) also gives a good estimate of how 
well one can do near fcmax in a scale dependent measurement. Upper right panel: Projected absolute error on Poe for the same 
cases, except the upper black (solid) curve is for perfectly observed velocity divergence (or two perfect tracers with different 
bias). Lower left (right) panel: Projected error on /3 (a), which has no cosmic variance limit. 
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on a. a can be measured very precisely, e.g., even for the simple split of the EUCLID-like sample in half, a will 
be measured to ~ 0.1% for fcmax — 0.05 ftMpc"^, and ~ 1% for fcmax ^ 0.01 ftMpc"^. This will measure the scale 
dependence of a and thus allow a very precise test of the onset of non-linear biasing. 

2. Higher density, smaller volume space mission 

Another option to consider is a space mission with a somewhat smaller volume, but higher density of galaxies. 
For example, one of the SNAP satellite designs would allow for a 10000 sq. deg. survey at 0.9 < z < 1.6, with 
density n = 0.004 ( ft.-^ Mpc) " and mean bias b — 2.3 (D. Schlegel, private communication). The higher density 
than EUCLID is what we look for to make the multiple-tracer method powerful. We show projections in Fig. [7] As 
an example, we split the SNAP sample into b — 2.8 and b = 1.8 subsamples. The improvement in Pgg measurement is 
substantial, a factor of ~ 2 — 5 in effective volume, depending on /cmax- The enhanced density of SNAP over EUCLID 
would roughly cancel the factor of 2 in greater sky area in the EUCLID survey (note, however, that the EUCLID 
survey still comes out ahead overall because it covers a broader redshift extent). 

3. SDSS-III/BOSS 

Fig. [5] shows the achievable gains for survey configurations motivated by the SDSS-III/BOSS survey (s^- BOSS 
will cover 10000 sq. deg. with central redshift z = 0.5 and maximum z ~ 0.7, giving comoving volume 6 cubic Gpc/h. 
The planned galaxy sample has b ~ 1.9 and number density n = 0.0003 { h~^ Mpc) . Not coincidentally, the form of 
Fig. [5] is very similar to Fig. [S] The only quantity that really matters for these calculations is the basic ratio of the 
noise level to power spectrum level, and the optimization of these surveys for BAO means that this ratio will be very 
similar. The BOSS S/N is somewhat lower than EUCLID, so BOSS is even farther from achieving the best possible 

measurement of Pgg. If we hypothetically split the sample into a 6 = 1.6, n = 0.0002 (ft^^Mpc) and b = 2.5, 

n = 0.0001 (ft^^Mpc) ^ subsamples (chosen for illustration, not because this is known to be possible), in order to 
exploit the method of ill, the gains are small ~ the noise is too large. The different panels of Fig. [8] show how well 
one can do in an absolute sense with a 6 cubic Gpc/h survey, on Pgg, P, and a. 

4. Present SPSS 

To put the future in context, it is interesting to see what we could do with current SDSS data. First, we take 
the LRGs, which we assume cover 2 cubic Gpc/h centered at z ~ 0.33 (10000 sq. deg.) with a density n = 
0.000068 ( ft- 1 Mpc) ~ , and b = 1.9. Fig. [9] shows our usual comparison between the Pgg measurement obtained 
from the LRGs and the ideal measurement. In this case, noise substantially degrades even the single-tracer version 
of the measurement. Splitting the sample into high and low bias parts (we assume b — 1.5 and 6 = 2.3) results in a 
completely negligible improvement, and adding the LRGs on top of a perfect unbiased tracer does not improve the 
results much (beyond the unbiased tracer results). Figure [9] also shows the absolute results for Pgg, (3, and a. While 
the Pgg and /3 figures are largely provided as a reference for how well a single-tracer can do on these quantities. The 
lower right panel of Fig. [S] shows that the a measurement for the sample split in half could be an interesting exercise. 
The bias ratio can be measured to about 2% at fc ~ 0.05 ft Mpc~^, or 1% at fc ~ 0.1 ftMpc""'^. This should be precise 
enough to place interesting constraints on the form of non-linear bias. 

The SDSS main galaxy sample is generally not as good for LSS as the LRGs, because it covers only low redshift 
and thus much less volume; however, it does have the higher number density that we look for to make the multi- 
tracer method powerful. [2J| provide a realistic split into blue and red samples with bias ratio 0.7, each with 
n ~ 0.0044 (ft~^ Mpc)~ at z < 0.1 (center of volume z ~ 0.08). We take 6 = 0.9 and b = 1.3. This sample covers 
only 0.026 cubic Gpc/h. Fig. [TO] shows that the multi-tracer approach to measuring Pgg is quite effective at this 
number density. Unfortunately, the upper right panel of Fig. [TO] shows that the absolute measurement still is not very 
good, requiring one to go to fcmax > 0.1 just to get a 3 — cr level of detection, and to fcmax > 0.3 to measure Pgg to 
10% (which is surely not possible at low z because of non-linearity). The situation with (3 is similar, but it is again 
interesting to note, in the lower right of Fig. [TO] that the bias ratio for the realistic split sample can be measured to 
a few percent. One should keep in mind however that at high fc the tracers are not maximally correlated and this 
weakens the constraints. Cosmological simulations will be needed to address this in more detail. 
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FIG. 7: Upper left panel: Projected error (inverse variance) on the normalization of Pgg = f^Pm, relative to the error one 
could achieve by observing velocity divergence directly (or using two perfect tracers with different bias), as a function of fcmax, 
for the SNAP mission-motivated example described in the text. The black (solid) line shows SNAP-like galaxies treated as a 
single tracer. The magenta (dotted) and blue (dot-short-dashed) lines show single perfectly sampled tracers with b = 2.3 and 
6=1, respectively. Green (short-dashed) shows the SNAP-like galaxies split into two groups with different bias, as described in 
the text. Red (long-dashed) shows the SNAP-like galaxies plus a perfect unbiased tracer, while cyan (dot-long-dashed) shows 
a density n = 0.03 {h~^ Mpc) plus an unbiased tracer. Where there are two lines of the same type, the upper one uses a 

genuinely perfectly sampled tracer, while the lower one shows density ft — 0.03 {h~^ Mpc) ^. The other panels are related to 
the upper left as in Fig. [B] 
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FIG. 8: Upper left panel: Projected error (inverse variance) on the normalization of Pgg = f^Pm, relative to the error one 
could achieve by observing velocity divergence directly, or with two perfectly sampled tracers, similar to Fig. |6] except this is 
for the SDSS-III/BOSS-galaxy-related example described in the text. The black (solid) line shows BOSS galaxies treated as 
a single tracer. The magenta (dotted) and blue (dot-short-dashed) lines show single perfect tracers with b = 1.9 and b — 1, 
respectively. Green (short-dashed) shows the BOSS galaxies split into two groups with different bias, as described in the text. 
Red (long-dashed) shows the SDSS-III galaxies plus a perfect unbiased tracer, while cyan (dot-long-dashed) shows ten times the 
BOSS galaxy density plus a perfect unbiased tracer. Where there are two lines of the same type, the upper one uses a genuinely 
perfectly sampled tracer, while the lower one shows density n = 0.003 (/i~^Mpc) ^ (long-dashed and dot-short-dashed) or 
n — 0.03 {h~^ Mpc) ^ (dot-long-dashed). The other panels are related to the upper left as in Fig. [6] 
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C. Beyond redshift-space distortions 



So far, primarily for the purpose of transparency of the presentation, the only information we have used from galaxy 
surveys is the BAO scale and Pgg{z,k = 0.1 /i Mpc~^), as measured from redshift-space distortions, ignoring other 
effects on the galaxy power spectrum. Generally, however, several other obvious effects should be included: The power 
spectrum cannot be observed in Mpc units, so the conversion to velocity (redshift) coordinates, Av = H{a) a /S.x\\ , 
and angular coordinates, A6 = D^(a)~^ a Axj_, must be accounted for. The uncertainty in the component of these 
conversions that takes one from angular coordinates to velocity coordinates, H{a)DA{a), is often called the Alcock- 



Paczyhski effect [17[. In addition to the transformation to observable coordinates, cosmological parameters generally 



affect the scale dependence of the power spectrum through the transfer function between early-time perturbations 
from inflation and the linear theory perturbations that we see, e.g., the classic F = rimfih scale for the turnover of 
the power spectrum. In Figure [11] we show the constraining power for the same general 3 / 4-sky survey up to z^ax 
that we have discussed before, except now computing the Fisher matrix using the full parameter dependence of the 
power spectrum. Note that the BAO measurement is automatically included within this framework, although in a 
somewhat weaker than usual form because of our relatively conservative fcmax cut. The extra parameter dependence 
leads to typically a factor 3-5 larger FoM than one would find using pPm and BAO constraints alone. There is a 
new feature at high S/N, however, in that we no longer find any cosmic variance limit on the constraining power! The 
measurement can be improved arbitrarily much if the S/N can be improved, but to achieve these improvements the 
number density of tracers must greatly exceed current designs. 

In Fig. [11] we also explore the complementarity of weak gravitational lensing measurements, as represented by the 
DETF Stage IV space (SNAP) or ground (LSST) projections (we use their optimistic versions). We see that these 
experiments, in combination with Planck (and all Stage II experiments, although these play a minor role), can achieve 
an FoM improvement of about a factor of 10 over the baseline. Redshift surveys are at best comparable if they are lim- 
ited to BAO, but they can potentially greatly exceed this if all information at fcmax('Z) ~ 0.1 [D (z) /D (0)]~^ /iMpc"'^ 
can be used. 



III. NON-GAUSSIANITY 



Many of the non-standard models of how the primordial structure was seeded predict non-Gaussianity of local 
type, <^ = (j) + fj^h4''^, where $ is the gravitational potential in the matter era and is the corresponding primordial 
Gaussian case. The effect of this type of initial conditions on the galaxy power spectrum was recently computed 
by [l^ and further investigate by [2^, HE [13, . The effect was used to place observational limits on /nl by [29[ . 
We can easily generalize the (l^ multi-tracer error on the non-Gaussianity parameter /nl to the more transparent 
approach of this paper. In that case, where redshift-space distortions have so far been ignored, the perturbations of 
galaxy type i can be written 

Sgi = {h + Cif^i^) 5 + ni (16) 

where Cj is an ideally known coefficient characterizing the galaxy type's response to non-Gaussianity (called A6(fc) 
by with A&i(fc) = 3(6; — pi) Sc^mHQ /c^k^T (k) D{z) from a model for clustering of halos - see [l4| for an 
explanation). The model we will use for perturbations of two types of galaxies is 

<5gi = 6 (h-ci/nl) (5-^711 , (17) 

where /nl = hh/b, and 

Sg2 ^b(a + C2/nl) S + n2 . (18) 

Note that the quantity we can measure precisely is not really /nl but instead /nl /nl/&- If one is trying to 
calculate the detection limit, as in [T^ . then the error on b is irrelevant, but for actual constraints on its value once it 
is detected the degeneracy between b and /nl becomes relevant. Like [l4|, we start by assuming that a is measured 
precisely from relatively small scales where the non-Gaussian effect is small, and Ci are known, so we only have two 
free parameters, P = 2 b^ ^^^1 /nl- The Fisher matrix for a single mode can be inverted to give the error on 
/nl, in the low noise, small /nl limit: 
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where, again, Xij = Nij /bibjPm- This result is the same as i n UM . once we account for differences in definitions. 

Estimates of the improvement factor for specific surveys in [14| are somewhat optimistic, because the improvement 
factor was evaluated using the signal-to-noise ratio a,t k = 0.01 /iMpc^^, near the peak of the power spectrum. This 
will give the correct improvement factor for a measurement using modes on that scale, but, as shown by p^, it is 
larger scale modes, where the signal-to-noise ratio for fixed noise will be lower, that will give the most interesting 
constraints on /nl- It is also not obvious when the assumption that P and a are perfectly measured using higher k 
data is safe. It seems useful to recalculate the constraints in the slightly more general approach we used for redshift- 
space distortions, i.e., including a as a free parameter and integrating over k up to some fcmax- The results are shown 
in Figs [12] and [131 Fig- [HI shows the same SDSS-III scenarios discussed for constraints on (3. For the single galaxy 
type case, we expect an error b a^^^ ~ 21, as long as we can use modes up to /Cmax ~ 0.05ft.Mpc~^. Note that the 
need to marginalize over the amplitude leads to a ~ 33% expansion in the errors. The improvement when the sample 
is split into high and low bias parts is modest, with b Uj^^ ~ 17. This is consistent with the finding of [l3|- Adding a 
perfect unbiased tracer helps substantially more, allowing a measurement to b o'^^^ ~ 7.3, again consistent with the 
improvement factor of |14| . The marginalization over nuisance parameters has led to a ~ 45% increase in the error. 

Fig [13] shows the /nl constraints for the EUCLID-like scenario, as discussed for /3. With only a single tracer, we 
expect an error ba^^^ ^ 2.0, as long as we can use modes up to fcmax ~ 0.05/iMpc~^. A simple split of the sample 
could improve this significantly, to 1.4 (note that [l^ assumed a factor of two higher galaxy density for his high 
redshift example, corresponding to the original SPACE plan, and so found a somewhat larger enhancement factor). 
The enhancement when a perfect unbiased tracer is added is dramatic, ba; ~ 0.54. In order to produce a solid 
detection for /nl < 1, we would need either a higher number density, e.g., like the SNAP sample we discussed above 
(except over a larger volume), or we would need to go to higher redshift where there is more volume, as discussed by 

MMM. 



IV. DISCUSSION AND CONCLUSIONS 



Galaxy redshift surveys allow one to measure the 3-dimensional distribution of galaxies, which contains an enormous 
amount of information compared to the 2-dimensional surveys such as cosmic microwave background anisotropics or 
weak lensing. Their main limitation is that we cannot directly determine the bias b which relates galaxies to the dark 
matter. In this paper we propose a new method to measure the redshift space distortion parameter (3 = f /b, where 
b is the tracer bias, / = din D/ din a, D is the growth factor and a the expansion factor. The method is based on 
angular dependence (with respect to the line of sight) of relative clustering amplitude of two tracers with different 
bias. In such a ratio the dependence on the mode amplitude vanishes and the method circumvents the usual cosmic 
variance limit caused by random nature of the modes. This allows one to determine the relative bias to an accuracy 
that is only limited by the number density of tracers. In principle this allows one to determine the velocity divergence 
power spectrum Pgg as precisely as one would by directly measuring the velocity divergence itself. This, in turn, 
would allow for a high precision determination of growth of structure as a function of redshift, as encoded in fD. A 
high density redshift survey of galaxies would allow for a more precise determination of dark energy parameters, as 
encoded in DETF Figure of Merit, than any other currently proposed Stage IV experiment. 

To achieve these gains one needs low noise, which in turn requires high density of galaxies for both samples. Roughly, 
the method requires nPg ^ 1, where n is the number density and Pg{k) is the galaxy power spectrum amplitude. 
However, most of the existing and planned BAO redshift surveys strive for nPg ^ 1, since one does not gain much 
in the galaxy power spectrum error once the field is significantly oversampled (i.e. fiPg ^ 1). Most redshift surveys 
only try to measure the shape of the power spectrum (in search of baryonic acoustic oscillations and other features) 
and not its amplitude. As long as survey speed is limited by the rate of observing individual objects (e.g., doubling 
the number density means halving the volume), the correct strategy will probably be to observe the full sky at the 
usual optimal BAO density, and then to go back to observe at higher density if possible. 

In addition to very precise measurements of /3, this multiple tracer method will allow precise consistency checks of 
the model. As one example, we can measure the bias ratio as a function of scale without cosmic variance and therefore 
study the onset of nonlinear bias. As another example, we can test the redshift space distortion model by fitting for 
a parameter /?' of the distortion model (l -I- -I- /3'/i^)j fo verify that /3' = 0. This new parameter (3' can also be 
measured to the same in principle bottomless level of precision as (3. 

In the short term, the technique in this paper will be most useful if methods are found that can achieve high 
effective number density without sacrificing volume, e.g., 21cm emission mapping (20j . Suppose one had a single 
perfectly sampled, high resolution field, say, from 21 cm emission mapping [20|. One can always create a second, 
biased field by some simple non-linear transformation, like squaring the original field, or taking just the high sigma 
peaks of the field. It may seem impossible that this could gain one any information, and that would be true if one 
was able to theoretically describe and exploit the field on all scales, including using higher order statistics; however. 
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when one is otherwise hmited to only very large scales this approach can add information, because it makes selective 
use of small-scale information, exploiting our understanding from renormalization of perturbation theory that local 
non-linear transformations of the density field can do nothing but change the bias parameters on large scales [l^ [l^ . 
Another way to see that this can work is simply to note that the two types of galaxy fields we have been discussing are 
simply two different non-linear transformations by nature of the same underlying density field, i.e., there is no reason 
we cannot do the same thing artificially. One would need to keep in mind in this approach that the transformation 
is applied to the redshift-space field, unlike standard bias which is a transformation of the real-space field. For the 
/nl application this is unlikely to be a fatal problem, but it may be for the redshift-space distortion measurement 
discussed in this paper. 

Our results suggest that there is much more to be gained by oversampling the density field and that one can 
determine the growth of structure as a function of redshift to a higher accuracy than previously believed. It is thus 
worth revisiting the planning of redshift surveys in light of these results. There is every reason to believe that these 
measurements can be made very robustly. We have boiled the constraint from a galaxy survey at a single redshift 
down to one number, but we will still have the full scale and angular dependence of the power spectrum, along 
with higher order statistics, to help us verify that there is nothing about the modeling that we do not understand. 
In addition, a realistic assessment would also need to address the issue of stochasticity between the tracers with 
numerical simulations, which would determine the value of kmax beyond which the two tracers are not strongly 
correlated anymore. Preliminary results suggest this is around kmax ^ 0.1 /iMpc~^ [s^l, but a more detailed analysis 
of multi-tracer method with simulations is needed, similar to the recent analysis of single tracer method [13]. We 
leave this subject to a future work. 

Another promising direction is to combine our no cosmic variance measurement of /3 = //6 with an analogous 
measurement of bias b using a comparison of weak lensing and galaxy clustering [3l|, to derive f{z) alone without 
any cosmic variance limitation. We leave these subjects to a future work. 

We than Ue-Li Pen, Nikhil Padmanabhan and Martin White for helpful discussions. U.S. is supported by the 
Packard Foundation, DOE and Swiss National Foundation under contract 200021-116696/1. 
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FIG. 9: Upper left panel: Projected error (inverse variance) on the normalization of Pgo = f'^Pm, relative to the error one could 
achieve by observing velocity divergence directly, or with two perfectly sampled tracers, similar to Fig. (6] except this is for the 
current SDSS LRGs. The black (solid) line shows LRGs as a single tracer. The magenta (dotted) and blue (dot-short-dashed) 
lines show single perfect tracers with b — 1.9 and b — 1, respectively. Green (short-dashed) shows the LRGs galaxies split into 
two groups with different bias, as described in the text. Red (long-dashed) shows the LRGs plus a perfect unbiased tracer. The 
other panels are related to the upper left as in Fig. [B] 
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FIG. 10: Upper left panel: Projected error (inverse variance) on the normalization of Pgg = f^Pm, relative to the error one 
could achieve by observing velocity divergence directly, or with two perfectly sampled tracers, similar to Fig. |6l except this 
is for the current SDSS main galaxies. The black (solid) line shows SDSS galaxies as a single tracer. The magenta (dotted) 
and blue (dot-short-dashed) lines show single perfect tracers with 6=1.1 and 6=1, respectively. Green (short-dashed) shows 
the main galaxies split into two groups with bias b — 1.3 and b — 0.9, both with n — 0.0044 {h~^ Mpc) , as described in the 
text. Red (long-dashed) shows the main galaxies plus a perfect unbiased tracer. Cyan (dot-long-dashed) shows two tracers 
with b = 1.3 and b = 0.9, both with n = 0.044 ( Mpc) ^. The other panels are related to the upper left as in Fig. |6l 
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FIG. 11: Improvement in DETF FoM when all of the information in the redshift-space power spectrum is used rather than 
just BAO and fPm constraints, for fcmax(z) = 0.1 [D (z) /D (0)]~^ /iMpc"^ (relative to Planck+Stage II, and for 30000 sq. 
deg., as usual). Green (dashed) hues show single S/N=l (at k = Q.4/iMpc~^) tracers with b = 1 (upper) and b — 2 (lower). 
The black curves show both types of galaxy together, each with (bottom to top) S/N=l, 3, 10, 30, 100, 1000 (the constraining 
power increases without bound as the S/N improves). Blue (dotted) shows BAO alone, with the lower curve using the above 
fcmax and the upper curve using smaller scales as described in [22 1 (the rest of the curves effectively use the weaker version). The 
red (long-dashed) curves include weak lensing constraints, with the lower two (horizontal) lines including no redshift survey 
(lowest is approximately as SNAP weak lensing, 2nd lowest is LSST weak lensing), while the upper two include the b = 2-only 
(weakest) redshift survey. 
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FIG. 12: Projected /nl errors (inverse variance) as a function of fcmax, for the SDSS-III/BOSS galaxy example described in 
the text. Strictly speaking, this is not (t/„l but instead b o-j:^^, where /nl = 6~^/nl, but the distinction is irrelevant for a 
rough detection limit. In each case, the lower curve of a type is marginalized over P and a, while the upper curve is not. Black 
(solid) assumes BOSS alone, treated as a single tracer. Green (short-dashed) curves show the case when the BOSS galaxies 
are split into high and low bias subsamples, as described in the text, while magenta (dotted) curves are the same except the 
transfer function is ignored when computing the non-Gaussian effect, restricting it to large scales. Red (long-dashed) show an 
infinite S/N, unbiased, tracer added to BOSS, while blue (dot-short-dashed) show the same with no transfer function effect 
in the calculation of the non-Gaussian signal. Cyan (dot-long-dashed) shows a perfect unbiased tracer added to a tracer with 
b = 1.9 and ten times the BOSS number density. The yellow (short-dash-long-dash) curves do not show errors on /nl, the 
upper curve shows {P/2N) /lO, where P is the power spectrum and is the noise for the full set of BOSS galaxies, and the 
lower curve /nlc (fc) /b where c(fc) is our best guess at the BOSS galaxies' response to non-Gaussianity, b is their bias, and 
/nl = 30. 
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FIG. 13: Projected /nl errors (inverse variance) as a function of femax, similar to Fig. 1121 for the EUCLID-like scenario. 
In each case, the lower curve of a type is marginalized over P and a, while the upper curve is not. Black (solid) assumes 
EUCLID alone, treated as a single tracer. Green (short-dashed) curves show the case when the EUCLID galaxies are split 
into high and low bias subsamples, as described in the text, while magenta (dotted) curves are the same except the transfer 
function is ignored when computing the non-Gaussian effect, restricting it to large scales. Red (long-dashed) show an infinite 
S/N, unbiased, tracer added to EUCLID, while blue (dot-short-dashed) show the same with no transfer function effect in the 
calculation of the non-Gaussian signal. Cyan (dot-long-dashed) shows a perfect unbiased tracer added to a tracer with b — 2 
and ten times the EUCLID number density. The yellow (short-dash- long-dash) curves do not show errors on /nl, the upper 
curve shows {P/2N), where P is the power spectrum and A*" is the noise for the full set of EUCLID galaxies, and the lower 
curve 1000/NLc(fc) /b where c{k) is our best guess at the EUCLID galaxies' response to non-Gaussianity, b is their bias, and 
/nl = 1. 



